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Gravity can become strong at the TeV scale in the theory of extra dimensions. An effective 
Lagrangian can be used to describe the gravitational interactions below a cut-off scale. In this 
work, we study the diphoton production in 77, pp, and e + e~ collisions in the model of low scale 
gravity. Since in the standard model photon-photon scattering only occurs via box diagrams, the 
cross section is highly suppressed. Thus, photon-photon scattering opens an interesting opportunity 
for studying the new gravity interaction, which allows tree-level photon couplings. In addition, we 
also examine the diphoton production at hadronic and e + e~ colliders. We derive the limits on the 
cut-off scale from the available diphoton data and also estimate the sensitivity reach in Run II at 
the Tevaton and at the future linear e + e~ colliders. 



I. INTRODUCTION 

Recent advances in string theories suggest that a special 11-dimension theory (dubbed as M theory) Jl| may be 
the theory of everything. Impacts of M theory on our present world can be studied with compactification of the 
11 dimensions down to our 3 + 1 dimensions. The path of compactification is, however, not unique. In this multi- 
dimensional world, the standard model particles live on a brane (3+1 dim) while there are other fields, like gravity 
and super Yang-Mill fields, live in the bulk. The scale at which the extra dimensions arc felt is unknown, anywhere 
from TeV to Planck scale. Recent studies Q show that if this scale is of order TeV and there are gauge and fermion 
fields living in the bulk that correspond to the Kaluza-Klein excitations of the gauge and fermion fields of the SM, 
early unification of gauge couplings can be realized below or even much below the original GUT scale. This is possible 
because the extra matters in the bulk accelerate the RGE running of the gauge couplings, which then change from 
logarithmic evolution to power evolution. Supersymmetry model building is also an active area in the framework of 
extra dimensions || . Apart from the above, radical ideas like TeV scale string theories were also proposed 0] . 

Inspired by string theories a simple but probably workable solution to the gauge hierarchy was recently proposed 
by Arkani-Hamed, Dimopoulos and Dvali (ADD) ||. They assumed the space is 4 + n dimensional, with the SM 
particles living on a brane. While the electromagnetic, strong, and weak forces are confined to this brane, gravity 
can propagate in the extra dimensions. To solve the gauge hierarchy problem they proposed the "new" Planck scale 
Ms is of the order of TeV in this picture with the extra dimensions of a very large size R. The usual Planck scale 
Mq = 1/ \/Gn ~ 1-22 x 10 19 GeV is related to this effective Planck scale Ms using the Gauss's law: 

R n M£+ 2 ~ Mq . (1) 

For n = 1 it gives a large value for R, which is already ruled out by gravitational experiments. On the other hand, 
n = 2 gives R < 1 mm, which is in the margin beyond the reach of present gravitational experiments. 

The graviton including its excitations in the extra dimensions can couple to the SM particles on the brane with 
an effective strength of 1/Mg (instead of 1/Mq) after summing the effect of all excitations collectively, and thus 
the gravitation interaction becomes comparable in strength to weak interaction at TeV scale. Hence, it can give 
rise to a number of phenomenological activities testable at existing and future colliders []fj|-p3|. So far, studies show 
that there are two categories of signals: direct and indirect. The indirect signal refers to exchanges of gravitons 
in the intermediate states, while direct refers to production or associated production of gravitons in the final state 
P,p| J^Jl^,pO|,pT|] . Indirect signals include fermion pair, gauge boson pair production, correction to precision variables, 
etc. [pp]9 Hl 7] , |19| , p2| , p3| . There are also other astrophysical and cosmological signatures and constraints [p3. 

Processes that only occur via loop diagrams in the SM are especially interesting if the low scale gravity allows tree- 
level interactions. In the SM, the lowest order photon-photon scattering can only take place via box diagrams of order 
a 2 (on amplitude level) (25) and, therefore, is highly suppressed. Thus, photon-photon scattering opens an interesting 
door for any tree- level photon interactions. Even if such new interactions are much weaker than the electroweak 
strength, these tree- level diagrams are only of order a new . It stands a good chance that these new interactions can 
beat the standard model. In the framework of ADD, photons can scatter via exchanges of spin-2 gravitons in s-, t-, 
and u-channels and the most important is that the coupling strength can be as large as the electroweak strength. In 
this work, we shall study the photon-photon scattering 77 — > 77 and demonstrate that it provides a unique channel 
to identify the low scale gravity interactions. Other interesting processes of the same category is 77 — > vv and the 
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cross-channel, 7^ — > 7^, both of which do not have any tree- level contributions in the SM |26|]. We shall not pursue 
these two further in this paper. 

Similarly, a pair of gluons can scatter into gluons or photons via exchanges of gravitons, the latter of which is our 
attention at hadron colliders. The lowest order gg — > 77 scattering occurs via a s-channel exchange of graviton in the 
low scale gravity model whereas it has to be via box diagrams in the SM. Thus, the new gluon scattering will give 
rise to anomalous diphoton production, in addition to the qq — ► G — > 77 channel, at hadron colliders. However, the 
tree-level SM qq — ► 77 presents a large irreducible background, not to mention the jet-fake background. This makes 
the diphoton production at hadron colliders not as attractive as in 77 and e + e~ colliders as a probe to the low scale 
gravity model. For completeness we also study the diphoton production at e + e~ colliders. 

The organization of the paper is as follows. In the next section, we compare the photon-photon scattering cross 
section between the SM and the low scale gravity. In Sec. Ill, we calculate diphoton production at the Tevatron and 
obtain the present limit on the cut-off scale Ms using the diphoton data, and then estimate the sensitivity reach at 
the Run II. In Sec. IV, we repeat the same exercise at e + e~ colliders and obtain the limits using the diphoton data 
from LEPII, and estimate the sensitivity reach at the future linear e + e~ colliders. We shall then conclude in Sec. V. 



II. PHOTON-PHOTON SCATTERING 



We concentrate on the spin-2 component of the Kaluza-Klein (KK) states, which are the excited modes of graviton 
in the extra dimensions. The spin-0 component has a coupling to the gauge boson proportional to the mass of the 
auge boson in the unitary gauge, which means it has a zero coupling to photons. We follow the convention in Ref. 
There are three contributing Feynman diagrams for the process 77 — > 77 in the s-, t-, and u-channels. The 
amplitudes for 7(^1)7(^2) - * 7(^1)7(^2) are given by: 



iMi = B^'^'ipx - ki) e p ( Pl ) e°{ P2 ) e a (k 1 ) e^{k 2 ) 

X[~Pl ■ k x C^.pa + ZV, pQ Ol,-fcl)] [~P2 ■ k 2 + Dp V ,a(3 (j>1 , ~k 2 )] 

iM 2 = iMi(ki <-> k 2 ) , 

^3 = ~\d{s) B^' v '{j> x + P2 ) eP{p x ) e°(p 2 ) e a (k x ) e {k 2 ) 

x[pi -p 2 Cpv^c + D f j, l/i p tT {p 1 ,p 2 )] [hi ■ k 2 C^ v ' >a p + Du'v^^-ku-ki)] , 



(2) 
(3) 



(4) 



where k = \Z16ttGn and Bp iV ^ pa {k), Cp V ^ pa and D tiUtPa -(pi,p 2 ) can be found in Ref. 



The propagator factor 



D(s)=J2 k i/(s- 



ie), where k sums over all KK levels. After some tedious algebra the square of the amplitude, 



summed over final and averaged over the initial helicities, is surprisingly simple: 



|D(s)| 2 (s 4 + t 4 + u 4 ) , 



(5) 



where we have taken M| ^> s, \u\ and in this case the propagator factor D(s) = D(\t\) = D(\u\) 0, which is given 
by 



K 2 \D(s)\ 



16tt 
Ml 



x T 



where the factor T is given by 



T =< 



log 
2 

n-2 



Mi 



The angular distribution is 



da(77 — * 77) TTi 



d\ cos t 



Ml 



for n = 2 
for n > 2 . 



l + i(l + 6cos 2 



(6) 



(7) 



(8) 



where | cos6*| is from to 1. Since the cross section scales as s 3 /M|, which implies larger cross sections at higher y^s. 

The SM background calculation is well known and we do not repeat the expressions here. We used the results in 
Ref. [^5| with the form factors from Ref. |27]] . The process is via box diagrams with all charged fermions and the W 
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boson in the loop. At the low energy, the fermion contribution dominates, but once -y/i gets above a hundred GeV 
the W contribution becomes more important and completely dominates at higher y/s. We show the cross sections in 
Fig. 0(a). This SM cross section decreases gradually when y/s is above 500 GeV. In contrast, the low scale gravity 
interactions give a monotonically increasing cross section. For n = 2 and Ms = 4 TeV the cross-over is at about 
■s/s = 600 GeV. We notice that the signal cross section does not decrease very rapidly with n, unlike the production 
of real gravitons P,[18[ . 

In Fig. |l](b), we show the angular distribution for the low scale gravity and for the SM. The signal has a relatively 
flat distribution, as can be easily deduced from Eq. (||). The ratio of the cross section at |cos#| = to that at 
| cos6>| = 1 is only 9/16. On the other hand, the SM background is very steep around | cos#| = 1, and that is why a 
cut of | cos0| < cos 30° is imposed to reduce the background. 

Monochromatic photon beam can be realized using the back-scatter laser technique |2£| by shining a laser beam onto 
an electron or positron beam. A linear e + e~ collider can be converted into an almost monochromatic photon-photon 
collider, with a center-of-mass energy about 0.8 of the parent e + e~ collider and with a luminosity the same order as 
the parent, i.e., as large as 50-100 fb _1 per year. Since the cross section for the SM is of the order of 10 fb, so there 
should be enough events for doing a counting experiment. A 5-10 % deviation from the SM prediction would be at a 
1.1 — 3.2<7 level. We use the 5% or 10% deviation from the SM as the criterion for sensitivity reach. The sensitivity 
reach at the 77 collider is shown in Fig. ||. The reach on Ms is about 5-8 (4.5-7.5) times of the center-of-mass energy 
of the collider for n = 2,4,6 using the 5% (10%) deviation criterion. As we shall see later, the sensitivity reach at 
photon-photon colliders is better than at e + e~ and much better than at hadron colliders. 



III. DIPHOTON PRODUCTION AT THE TEVATRON 



Diphoton production has been an interesting subject for CDF and DO. It can provide constraints on the (7977 
type contact interactions, anomalous 777 and 277 couplings. In the context of the low scale gravity, diphotons can 
be produced via quark- antiquark and gluon-gluon annihilation into virtual gravitons and the associated KK states. 
The gluon-gluon annihilation is very similar to the photon-photon scattering described in the last section. The main 
background is the SM lowest order process: qq — ► 77. [] 

There are two contributing subprocesses: 



dojqq -> 77) 
dcosO* 



1 



967TS 
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27^(2^(1 + 003^*)^ 



Ml 



(i-cos 4 r)^ 2 



(9) 



^(55^77) * S 3 2 4 2 

dcosfl* =5T2Mf (1 + 6C ° S d +C ° S 6)T ' (10) 

where the factor T is given in Eq. (0), and the 0* is the scattering angle in the center-of-mass frame and cos 9* is 
from —1 to 1. In qq — > 77, the effect of graviton exchanges first occurs in the interference term, which only scales as 
s 2 /Mg, and potentially more important than the square term of s 4 /M§ at s <C M§. 

Both CDF and DO |29| have preliminary data on diphoton production. We are going to use their data to constrain 
Ms- CDF has measured the invariant mass M 77 spectrum in the region 50 GeV < Af 77 < 350 GeV. However, since 
the data is only preliminary and in graphical form only, we can only use the reported number of events in the region 
M 77 > 150 GeV: 5 events are observed where 4.5 ± 0.6 are expected with an integrated luminosity of 100 pb _1 . 
This data, though without binning information, is sufficient to place a constraint on Ms, because the signal for the 
low-scale gravity does not appear as a peak in the M 77 spectrum but, instead, a gradual enhancement from about 
M 77 w 150 GeV towards higher M 77 . We use the Poisson statistics to calculate the 95% CL upper limit to the 
number of signal events V95, n using 



1 Since the lowest order diphoton production qq — > 77 is much larger than the box process: gg — > 77, we shall neglect the 
latter in considering the SM background. 
2 The number of signal events is N95 or less with 95% confidence. 
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p -(n B +N 95 ) y^nobs (wa+jVos)" 

0.95=l-e=l ^^2-5 — _»! (11) 

where jib = 4.5 is the expected number of background events and n b s = 5 is the number of observed events. We 
obtain iVgg = 6.61. We then normalized our calculation to the expected number of events (=4.5) after imposing the 
same selection cuts as CDF. With this normalization we can then calculate Ms, which gives a signal of 6.61 events in 
excess of the SM prediction. We obtain the 95% CL lower limit on M$: 

Tevatron Run I: M s > 0.91 TeV for n = 2 and 

M s > 0.87 TeV for n = 4 . 

For DO, however, the highest bin in the measured M 77 spectrum is 80-112 GeV. At such a low value, it is difficult to 
see the effect of (s 2 /M|). Thus, we expect the limit that would be obtained from the DO data is somewhat smaller 
than using the CDF data. 

Next, we estimate the sensitivity reach at the Run II of the Tevatron, assuming a luminosity of 2 fb _1 . The effect 
of the low scale gravity on the M 77 spectrum is shown in Fig. |^. It is easy to understand why the enhancement 
is more likely at the large M 77 . To estimate the sensitivity we divide the M 77 spectrum into bins: a bin width of 
100 GeV for bins in 200 GeV < M 77 < 500 GeV, and for 500 GeV to 1000 GeV we combine it into one bin only. 
This is to make sure that each bin should have at least a few events in the SM: see the first row of Table Q (we also 
use a selection efficiency of 50%.) For each bin we assume the SM prediction as the number of events that would be 
observed: n obs , and we calculate the number of events predicted by a Ms and n: n th . We then calculate the % 2 for 
this bin and sum over all bins, using 



X 2 (M s ,n)= £ \ 2 (n?-n?»)+2n?°]n(^) 



obs 

(12) 



i— bins 

The x 2 then gives a goodness of the fit for the value of Ms and n. The larger the \ 2 the smaller the probability that 
the corresponding value of Ms and n is a true representation for the data. To place a 95%CL lower limit on Ms a 
X 2 = 9.49 is needed for 4 degrees of freedom. The number of events in each bin for n = 2 and Ms = 1.5 — 2 TeV, and 
for n = 4 and Ms — 1.4 — 2 TeV with the corresponding x 2 are shown in Table ffl. We obtain a limit of 

Tevatron Run II: M s > 1.72 TeV for n = 2 and 

M s > 1.43 TeV for n = 4 . 

We verified that the binning is not important for the limit. We repeat the procedures using only one large bin from 
200 to 1000 GeV, and the 95% CL lower limit on M s becomes 1.73 (1.38) TeV for n = 2 (4). 

TABLE I. The number of events that would be observed in each bin of M 77 for the SM and for the low scale gravity at 
the Tevatron with <Js = 2 TeV and a luminosity of 2 fb _1 . The x 2 ls calculated assuming the SM prediction is what would be 
observed. The cuts imposed are: |jj 7 | < 1, pt 7 > 20 GeV, and a selection efficiency of 0.5 is assumed. 



model 


200-300 GeV 


300-400 GeV 


bin 

400-500 GeV 


500-1000 GeV 


x 2 


SM 


47.68 


11.98 


3.65 


1.81 




n = 2 
Ms = 2.0 TeV 
M s = 1.9 TeV 
M s = 1.8 TeV 
M s = 1.75 TeV 
M s = 1.7 TeV 
M s = 1.6 TeV 
M s = 1.5 TeV 


50.27 
50.82 
51.53 
51.96 
52.47 
53.75 
55.49 


14.40 
14.92 
15.59 
15.99 
16.45 
17.61 
19.25 


5.53 
5.93 
6.46 
6.78 
7.15 
8.09 
9.38 


4.84 
5.54 
6.45 
7.01 
7.66 
9.30 
11.58 


3.80 
5.24 
7.33 
8.70 
10.37 
14.87 
21.72 


n = 4 
Ms = 2.0 TeV 
M s = 1.9 TeV 
M s = 1.8 TeV 
M s = 1.7 TeV 
M s = 1.6 TeV 


48.24 
48.38 
48.54 
48.76 
49.07 


12.62 
12.77 
12.97 
13.24 
13.62 


4.22 
4.37 
4.55 
4.80 
5.15 


2.96 
3.28 
3.72 
4.38 
5.35 


0.64 
0.97 
1.49 
2.39 
3.89 
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M s = 1.5 TeV 
M s = 1.4 TeV 



49.52 
50.20 



14.16 
14.94 



5.65 
6.40 



6.87 
9.35 



6.53 
11.30 
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IV. DIPHOTON PRODUCTION AT E + E~ COLLIDERS 



We can use Eq. (M) with Q q = — 1 and multiply it by 3 to derive the expression for 



77: 



da(e + e — > 77) 2n 
dz s 



, 1 + z 2 
1-z 2 



4 Mi. { 



1 



(13) 



where z = | cos(9| is the polar angle of the outgoing photon and z ranges from to 1. 

The four LEP collaborations have been measuring the diphoton production e + e - — * 77 ||(| and using the data to 
constrain the deviation from QED and generic types of contact interactions of order 1/A n , n = 6,7,8. Since these 
contact interaction parameters 1/A n can be converted from the QED cutoff parameter A±, we shall stick with the 
QED cutoff parameter in the following discussion. The possible deviation from QED is usually characterized by a 
cutoff parameter A± corresponding to a modified angular distribution: 



da 

dz 



2-wa 2 1 + z 2 



1 



1 ± 



2K% 



(1-z 2 ) 



(14) 



where z = | cos^l and ranges from to 1. 

Each collaboration measured the cos 9 distribution and obtained the 95% CL limit on A± by varying r\ = 1 / A ± and 
maximizing the likelihood function. Since each experiment has their own procedures, we adopt a simple approach 
that takes their limits on A± and converts them into limits on Mg. Note that in Eq. (|l3|) the third term is suppressed 
relative to the second term, we can, therefore, just take the first and the second term, then it will look like Eq. (|TT 
The QED cutoff parameter A + is related to Mg by 



Ml 
T 



K 

2a 



(15) 



The limits from each LEP experiment and the corresponding limits on Ms are tabulated in Table [fi|. Note that 
we used only A + to calculate Mg. The limits on Ms is at most about 1.4 TeV for n — 2 and about 1 TeV for 
n — 4. The result for n — 2 is enhanced because of the logarithmic factor in T . Using the value of Ms ~ 1 TeV 
we can verify the ratio of the third term to the second term in Eq. ( |l3| ) and the third term is only about 2% of the 
second term. It justifies the approximation that we take only the first two terms of Eq. ([l3]). So far, the treatment 
is rather simple. A better limit can be obtained by combining the data on 77 = 1/A* from each LEP experiment. 
However, since some of the data on r\ are not given in detail, we can only combine those with a central value and an 
error. We have the following available: (i) OPAL (183 GeV): rj = (1.04 ± 1.34) x 1CT 10 GeV -4 , (ii) L3 (183 GeV): 

r] = (-0.59 ±i:i§) x 10 -10 GeV -4 , (iii) L3 (161,172 GeV): rj = (-0.77 -iil) x 10 -10 GeV -4 , and (iv) DELPHI 
(183 GeV): r; = (—1-4 ± 1.5) x 10~ 10 GeV - . We combine these data and assuming they are all gaussian we obtain 

77 = (—0.31 -0.73) x 10 -10 GeV -4 , the error of which is given in ler. From this i] the corresponding 95% CL limit on 
A± arc A + > 298 GeV and A_ > 279 GeV. We can see that the combined limit on A + is still not as good as the 
single limit from ALPEH (189 GeV) or OPAL (189 GeV). Once the data from each LEP experiment are given we 
can certainly improve the limit by combining them. Thus, for the present moment the best limit is from OPAL (189 
GeV): A+ > 345 GeV, which converts to M s > 1.38 (0.98) TeV for n = 2 (4). 

The behavior of the new gravity interactions at higher y/s can be easily deduced from Eq. (|l3|). The new interaction 
gives rise to terms proportional to s 2 /Mg and s 4 /Mg, which get substantial enhancement at large y/s: see Fig. f|(a). 
The angular distribution also becomes flatter because in the SM the distribution scales as (1 + z 2 )/(\ — z 2 ) whereas 
the terms arising from the new gravity interactions scale as (1 + z 2 ) and (1 — z 4 ), respectively, as shown in Fig. ||(b). 

Here we also attempt to estimate the sensitivity reach on the cut-off scale Mg at the future linear e + e - colliders. 
Since the cross section is of the order of 0.1 to 1 pb for yfs — 0.5 — 2 TeV, it corresponds to about 10 3 — 10 events 
for a mere yearly luminosity of 10 fb -1 . Thus, a 5% (10%) deviation from the SM prediction corresponds to a level of 
1.6cj — 5cr (3.2<7 — 10ct). In Fig. |s|, we show the sensitivity reach on Mg by requiring a 5% or 10% deviation from the 
SM prediction. The reach on Ms is about 3.5-5.5 (3-4.5) times of the y/s of the collider for the 5% (10%) criterion. 

TABLE II. The 95% CL limits on the QED cutoff parameter A± from LEP experiments Q and the corresponding 95% 
limits on Ms obtained using Eq. (|l5|). We show only the result of the highest energy of each experiment whichever available. 



95%CL limit on A+ and A_ 


95% CL limit on M s (TeV) 
n = 2 n = 4 


OPAL (Vs= 189 GeV): 


A+ > 345 GeV 


1.38 0.98 
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A_ > 278 GeV 




Pi F T PHT (+n — 1 K\ HpVI • 


A_ > 225 GeV 


n Q7 




L3 (y/s = 183 GeV): 


A+ > 262 GeV 
A_ > 245 GeV 


1.01 


0.74 


ALEPH (Vs = 189 GeV): 


A+ > 332 GeV 
A_ > 265 GeV 


1.32 


0.94 
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V. CONCLUSIONS 



Diphoton production at 77, pp, and e + e~ colliders provides useful channels to search for the presence of the low 
scale gravity interactions, which are the effects of allowing gravity to propagate in the extra dimensions. Photon- 
photon colliders are able to give the best sensitivity reach on the cut-off scale Ms of the low scale gravity model 
among the three. This is because 77 — > 77 can only occur via box diagrams in the SM while in e + e~ and pp collisions 
the tree-level contributions from the SM dominates. In addition to the total cross section, the angular distribution 
also serves as a tool to distinguish between the SM and the new gravity interactions, as seen in Fig. [l](b) and Fig. 
|(b). 

The present limit from the LEPII diphoton data is about Ms > 1.4 (1) TeV for n — 2 (4), and it is only Ms > 0.9 
TeV from the CDF diphoton M 77 data. The sensitivity reach in 77 collisions is about 5-8 times of ^/s 77 while it is 
only 3.5-5.5 times of the y/s at e + e~ collisions. At the Run II of the Tevatron, the reach is only about 1.7 (1.4) TeV 
for n = 2 (4). 

Finally, we emphasize the diphoton production at photon-photon colliders could provide a unique probe to the 
collider signature for the model of low scale gravity. 
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FIG. 1. (a) The total cross sections and (b) the differential distribution da/d\ cos# 7 | for 77 — > 77 for the low scale gravity 
model and for the SM. A cut of j cos # 7 j < cos 30° is imposed. 
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FIG. 2. The sensitivity reach on Ms versus ^/s 77 using the process 77 — > 77, by requiring the signal to be 5% or 10% of 
the SM prediction. A cut of | cos # 7 1 < cos 30° is imposed. 
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FIG. 3. The differential distribution d<r/c!M 77 versus M 77 for diphoton production at the 2 TeV Tevatron for the SM and 
for the low scale gravity with Ms ~ 1.5, 2 TeV and n = 2,4. Cuts of |rj 7 | < 1 and Pt 7 > 20 GeV are imposed. 
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FIG. 4. (a) The total cross section versus yfs and (b) the differential distribution da/d\ cos#| for e + e — > 77 for the SM 
and for the SM plus the new gravity interactions with n = 2, 4, 6 and Ms as shown. The j cos 8\ < 0.95 is imposed. 



13 



12 r- 




I — i 1 , 1 , 1 , 1 

500 1000 1500 2000 

Vs (GeV) 



FIG. 5. The sensitivity reach on Ms versus y/s using the process e + e — > 77, by requiring a 5% or 10% change from the 
SM prediction. A cut | cos 9\ < 0.95 is imposed. 
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